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We have experimentally observed the dynamics of an antiferromagnetic sodium Bose-Einstein 
condensate (BEC) quenched through a quantum phase transition. Using an off-resonant microwave 
field coupling the F = 1 and F = 2 atomic hyperfine levels, we rapidly switched the quadratic 
energy shift q from positive to negative values. At ^ = the system undergoes a transition from a 
polar to antiferromagnetic phase. We measured the dynamical evolution of the population in the 
F — l^rriF — ^ state in the vicinity of this transition point and observed a mixed state of all 3 
hyperfine components for ^ < 0. We also observed the coarsening dynamics of the instability for 
^ < 0, as it nucleated small domains that grew to the axial size of the cloud. 

PACS numbers: 03.75.Mn,67.85.De,67.85.Fg,67.85.Hj 



A quantum phase transition describes a many-body 
system whose ground state can be tuned through a point 
of non-analyticity [1]. Quantum gases afford the possi- 
bility to realize such phase transitions in the laboratory, 
as well as to explore the dynamical evolution of the state 
of the system by directly controlling the tuning parame- 
ters. In particular, spinor Bose-Einstein condensates pos- 
sess a vector order parameter with additional degrees of 
freedom relevant to this problem. By changing external 
fields dynamically, one can observe and quantify a host of 
non-equilibrium phenomena, including spin domain for- 
mation and aggregation, topological defect creation, and 
possible dynamic scaling laws [2]. 

In this work we examine a first order phase transi- 
tion associated with the quadratic energy shift, which is 
an essential parameter in spinor physics. For a spin-1 
BEC the spin-dependent Hamiltonian can be written in 
a mean-field representation as 
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where F,^^ are the vector spin-1 operator and its z- 
projection, respectively, n is the particle density, C2 the 
spin-dependent interaction coefficient [29], and q is the 
energy difference — ^o, where Ei is the en- 

ergy of the atomic level for the spin mp = i component 
of F = 1. The spin-dependent interaction coefficient C2 
arises from spin-changing collisions that can convert two 
mp = atoms into an mp = ±1 pair and vice- versa, a 
process constrained by the conservation of angular mo- 
mentum. It determines the nature of the ground state- 
antiferromagnetic for C2 > or ferromagnetic for C2 < 0. 
The quadratic energy shift q is usually determined by 
an external magnetic field B through the second-order 
Zeeman effect, and is cx B^. However, it can also be 
tuned using a microwave dressing field d, : a feature we 
exploit in this work to uncover a previously unexplored 
phase transition in an antiferromagnetic ^^Na spinor gas. 
While in general the levels shift independently, spin con- 
servation leads to the cancellation of the linear energy 



shifts, such that only the quadratic energy shift q is im- 
portant for spinor BECs. Therefore q plays the role of 
an external parameter and the combination of C2 and q 
realizes a rich phase diagram of possibilities [5]. Various 
quantum phases and dynamics have been observed for 
both C2 > and C2 < 

For an antiferromagnetic spinor BEC constrained to 
have zero net magnetization, the ground state solution 
is a nematic order parameter ^. It varies smoothly with 
q for all values except <7 = 0, which divides the phase 
diagram into two regions. For q > the ground state is 
a polar condensate consisting of a single component-the 
mp = spin projection that minimizes (F^). For g < 
the ground state maximizes the same quantity through a 
superposition of two components mp = ±1, a so-called 
antiferromagnetic phase [2] . The symmetry properties of 
the ground state therefore change discontinuously, defin- 
ing a first order phase boundary. 

Exactly at the phase transition point the many-body 
ground state is a condensate of boson pairs forming a 
spin singlet state, and possessing super-Poissonian spin 
fluctuations. Near this boundary the mean- field wave- 
function ^ undergoes collapse and revival dynamics trig- 
gered by quantum fluctuations [isl-llij. Controllably ac- 
cessing q = would restore the full S'^ symmetry of the 
nematic order parameter, which has unusual topological 
defects such as half-quantum vortices 16-18|. However, 
experimentally this requires low magnetic fields where 
spinor condensates are susceptible to ambient magnetic 
field noise that can mask interaction-related phenomena. 
In this letter, we investigated the dynamical instability 
ofanF = l^mp = antiferromagnetic sodium BEC 
that is rapidly quenched across the boundary from q > 
to q < 0. The quadratic energy shift q was tuned by 
an additional microwave dressing field that allowed us to 
access the q < regime. 

Microwave dressing of ferromagnetic ^^Rb has been 
used to tune spin mixing dynamics into the resonant 
regime in optical lattices [sj , and for the study of sponta- 
neous magnetization Unlike ferromagnetic systems. 
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in antiferromagnetic condensates one generally expects 
unmagnetized domains to form as a result of = ±1 
pair production. Pair formation dynamics have also been 
studied in antiferromagnetic F = 2 ^^Rb spinors where 
the sign of both C2 and q are reversed with respect to 
the F = 1 manifold [III, 0, [2o|, although the spinor 
lifetime is limited by hyperfine state changing collisions. 
The F = 1 spin system, by contrast, is intrinsically sta- 
ble and amenable to studies of long time scale dynamics 
in the transition region-in the present work we observed 
evolution times ranging from 30 ms to over 2 seconds. 

We prepared BECs with up to 3 x 10^ sodium atoms 
with a peak atom density no = 5.4 x 10^"^ cm~^ in an op- 
tical dipole trap created by a single focused far-detuned 
1064 nm laser beam. The measured axial trapping fre- 
quency was 8 Hz, with inferred radial frequencies of about 
600 Hz, which correspond to Thomas- Fermi radii of 4 jiia 
and 270 /im, respectively. The condensate is initially pre- 
pared in a magnetic trap from which it was transferred 
into the optical trap as described in earlier work [2L] . In 
order to create a pure = condensate, we adiabat- 
ically swept the frequency of an rf magnetic field across 
the rriF = —1 mp = transition at a bias magnetic 
field of 13 Gauss [22|. The bias was then ramped to a 
final value of 5 = 97 mG in 30 ms, which was defined as 
the starting point (t = 0) of our experiment. 

At this magnetic field, the quadratic Zeeman shift 
is qb = hx 2.5 Hz far smaller than the spin- 
dependent interaction energy C2no = /ixl30 Hz |30| . 
The radial Thomas-Fermi radius of our cloud, 4/im, is 
only a factor of 3 larger than the spin-healing length 
^5 = %/ \/2Mc2n 1.3 /im, which is the width of a typ- 
ical spin domain wall [23|. Therefore our experiment is 
performed mostly in a quasi- ID geometry, in which only 
longitudinal spin structures are expected. 

At t = we instantly turned on an oscillating mi- 
crowave field within tens of /is, far shorter than any dy- 
namical timescale relevant to the problem. We used mi- 
crowave powers between and 7.5 Watts, which allowed 
us to tune q from its initial value of to a final value 
Q = Qb Qm = -18.5 Hz, where qm is the quadratic 
energy shift due to the AC Stark shift caused by the 
microwave field [3]. This field created an instability in 
the initial tuf = spin state. Holding the condensate 
at a fixed value of q and varying the hold time after the 
quench was initiated, we could observe the temporal evo- 
lution of the fractional population in = 0. This cor- 
responded to a measurement of the squared amplitude of 
the z-component of the nematic order parameter. The 
relative populations in the tuf = — 1,0,+1 states were 
determined by Stern- Gerlach time-of-fiight images. After 
an expansion of 3-3.5 ms, we pulsed on a Stern-Gerlach 
field gradient for a duration of 2-4.5 ms perpendicular to 
the axial direction separating the 3 components spatially 
0. After a total time-of-fiight of 25 ms, we optically 
pumped the atoms into the F — 2 state and imaged 
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FIG. 1: Quenching dynamics. For different final values of ^, 
the rriF — ^ fraction (circles) is plotted versus time after the 
quench. The solid lines are fits of the mp — ^ fraction to a 
Sigmoid function. In a) no microwave dressing field is applied. 
This data shows the relaxation in the absence of a quench 
(g = /i x 2.5 Hz). The equilibration to a pure = ±1 cloud 
occurs after 2.1 s. In b) the gas is quenched to ^ = hx—3.2 Hz, 
showing that the population decay is faster. In c) the gas is 
quenched to q = h x —17 A Hz and rapidly reaches a quasi- 
equilibrium state containing all 3 spin components in roughly 
equal proportions, i.e. an order parameter delocalized over 
the S^ sphere. Each data point in the figure corresponded to 
a separate run of the experiment. 

them on the F = 2 ^ F' = 3 cycling transition ensur- 
ing an equal imaging sensitivity to each spin component. 
Examples of the temporal evolution for q = h x —3.2 Hz 
and q = h X —17AB.Z are shown in Figures [TJd) and c), 
respectively. They indicate that a pure m = conden- 
sate was unstable, evolving into a superposition of all 3 
components that preserved the overall zero net magne- 
tization. In order to quantify the g-dependence of the 
instability, we fitted data such as shown in Figure [T] to 
a Sigmoid function and determined both the cross-over 
time Ti/2 and the final saturation value /o,min where 
/o(^i/2) = 1/2(1 + /o,min)- We defined the measured 
instability rate to be r{q) = I/T1/2. 

At zero microwave power, the quadratic energy shift 
was <7b/(c2^o) = 0.02. Under these conditions, the gas 
was relatively stable against spin relaxation, i.e. the cre- 
ation of = ±1 pairs. This stability is a characteristic 
of quantum antiferromagnetism Q. For a homogeneous 
system and < q < C2n the uif = state is stabilized 
against the creation of magnon excitations (spin waves) 
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by an energy gap A ^/{e^ 

6fe = ^ 0. The k 
mode and for q 



Q){^k + <7 + 2c2n), where 
mode is the most unstable 
C2n, the energy gap 



<< C2n, the energy gap ^ ^/2qc2n = 
hx25 Hz for our parameters. 

As q ^ 0, quantum fluctuations destabilize the pure 
rriF = state: the fraction of atoms in the = ±1 
states should reach of order 1 within a time ~ 1.5 sec- 
onds p^]. This time scale is consistent with the slow 
rate of relaxation (~ 0.5 s~^) that we observed in our 
experiment (see Figure [T] a)). However, we cannot defini- 
tively rule out other mechanisms, including thermal fluc- 
tuations (our BEC had a thermal fraction of 40%) and 
imperfect transfer to the tuf = state. For hold times 
longer than 2 seconds the cloud separated into two non- 
overlapping rriF = ±1 spin domains along the long axis 
of our trap, an indication that small residual magnetic 
field gradients might have been present. 

For —C2n < q < the gas becomes unstable against 
pair formation due to the presence of an imaginary fre- 
quency F = Xm{A). As \q\ increases, the gas progres- 
sively evolves into a mixed state. In all cases the tempo- 
ral dependence of the tuf = fraction, /o(t), followed a 
backward S'-shaped curve that saturated at a value fo^min 
depending on q. The data sets could be roughly divided 
into 2 regions: for positive (Figure [T^)) and slightly neg- 
ative q (Figure [TJd)) we observed a slow decay to a value 
fo,min nearly zero. For this data, an examination of the 
spatial distribution of the three hyperfine components 
showed that the ±1 states had separated from one an- 
other along the axial direction. This is most likely due 
to the residual magnetic field gradients mentioned earlier. 

For more negative q (Figure [it)), however, the behav- 
ior was dramatically different-/o approached a final value 
of /o,min ^ 0.3 within a time as little as 30 ms and re- 
mained roughly constant over 200 ms [31]. Thus for the 
range of q explored in this work, the instability was ob- 
served to create a mixed state of all 3 components which 
appeared to be metastable on a timescale much longer 
than r~^. When accounting for residual thermal atoms, 
the estimated condensate fraction tuf = is slightly less 
than 0.3. 
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FIG. 2: Final rriF — fraction versus instability rate P. The 
open circles indicate data for which the final ttif = distri- 
bution consisted purely of thermal atoms. 



Figure [2] shows the final saturation value /o,min(^) plot- 
ted versus T{q). These data provide further evidence for 
the two regimes mentioned earlier. Above a critical in- 
stability rate of :^ 3 s~^, the final fraction /o,min was 
between 0.2 and 0.3, and the ttif = state maintained a 
significant presence in the cloud. With the exception of 
one data point, only for the very lowest instability rates 
< 3 s~^ was the steady state Bose-condensed spin dis- 
tribution consistent with a pure = ±1 spin mixture 
(open circles). 

We can qualitatively understand the instability in 
terms of a quantum rotor model that is valid in the single 
mode approximation [TBj . Such a model cannot describe 
spin domain formation, and therefore is inapplicable to 
the data for small \q\ discussed above, for which we ob- 
served spin segregation. However, it may provide useful 
insights into the short timescale behavior following the 
quench for larger \q\^ where a mixture of all 3 compo- 
nents was observed. The system is described by a single, 
macroscopic quantum rotor with angular momentum L, 
moment of inertia / = NTi^ / {c2n), and a potential energy 
term V ~ qNsiii?{0) for our parameters. ^ = describes 
a pure m = state, while < ^ < tt corresponds to the 
inclusion of ±1 pairs into the wavefunction. The sud- 
den quench transition from q > io q < ^ causes V to 
change sign, resulting in rapid dispersion of a wavepacket 
initially highly localized in angle near ^ = 0. The result- 
ing wavepacket dynamics are mostly classical in character 
and consist of a rapid dispersion in followed by sparse, 
periodic revivals at time trev [Hj . For our parameters 
trev ^ 350 seconds, considerably longer than our obser- 
vation time [15]. For short times, as in Figure [It), we 
observed only the rapid dispersion phase, and we inter- 
pret the measured value of ~ 0.3 for each spin component 
to be the result of wavepacket dispersion that tends to 
equalize the spin populations. 

The measured instability rate F has been plotted ver- 
sus the final quadratic energy q = qB-\-qM '^^ Figure [3] for 
q ranging from + 2.5 Hz to -18.5 Hz. The data shows a 
steep rise in the pair formation rate by a factor of nearly 
100 as — ^ varied from +0.02 to —0.15, indicating that 

C2no ' ^ 

we had crossed a phase boundary in the dynamical evo- 
lution of the system. Figure [3] also shows the predicted 
maximal instability rate for — C2no < <7 < from Bogoli- 
ubov theory for a uniform ttif = gas with the same 
average density (n) , Tuni f = \/\q\{q + 2c2(n)). This cor- 
responds to the formation of correlated pairs of atoms in 
a spatial mode with wavevector /c = 0, i.e. a homoge- 
neous rate of pair formation throughout the cloud [2J|. 
The homogeneous theory is in considerable disagreement 
with our data, which could be attributed to the finite 
size and 1-D geometry of our trap. This departure is 
consistent with earlier work on pair formation dynamics 
in F = 2 spinor condensates which highlighted resonant 
structures and the importance of the inhomogeneous den- 
sity profile in determining the modes that were populated 
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FIG. 3: Quenching through the quantum phase transition. 
The formation rate of = ±1 atom pairs is plotted versus 
the quadratic energy shift q (circles) determined by fitting 
the temporal evolution of the ttif — fraction to a Sigmoid 
function. The error bars give the statistical uncertainty in 
the fit. The instability rate dramatically increases below the 
transition point at ^ = 0. Also shown (dotted line) is the pre- 
dicted instability rate from Bogoliubov theory for a uniform 
gas. Inset shows the same data plotted on a semilog scale. 



Ill |26| . No clear indication of resonances were visible in 
our data. Moreover, this theory does not account for 
possible spin exchange processes between the condensate 
and the residual thermal cloud, which could play a role 
in our observations 
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In order to gain an intuitive understanding of the 
spatial dynamics one can think of the inhomogeneous 
mi? = condensate as a locally varying gain medium 
for the pair formation instability. For very short times 
after the quench, depletion of the gain can be neglected, 
and one may write the growth rate of mp = ±1 atom 
pair number for g < using a local density approxima- 
tion as V local = \/\q\{q + 2c2n(r))//i, where n(r) is the 
spatial density profile of the mp = cloud with the max- 
imum gain occurring at the cloud center. In this regime, 
the inhomogeneous gain acts as a nonlinear spatial mode 
coupler that converts energy from long to short wave- 
lengths, i.e. exhibits a tendency to nucleate small sized 
domains. These dynamics are not captured in the homo- 
geneous theory, but appear in our data. For data sets 
with q < hx -7 Hz, the mi? = ±1 atom distribution ap- 
peared initially as one or more small axial domains near 
the cloud center. As an example. Figure |4^)-e) show 
Stern- Gerlach images at various times after the quench 
for q = h X —17 A Hz. These domains appeared to co- 
alesce into a larger domain that grew in size with time 
until it became comparable to the axial Thomas-Fermi 
radius (see Figure Hf)). 

Once a substantial number of atom pairs have been 
created the mi? = condensate can be locally depleted 
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FIG. 4: Spatial dynamics of the instability. Absorption im- 
ages of the condensate taken at a time of flight (TOF) of 25 ms 
for a quench to ^ = x — 17.4 Hz for different hold times: a) 
15 ms, b) 20 ms, c) 25 ms, d) 30 ms, and e) 150 ms. From 
top to bottom the images show the mp — —1, mp — 0, and 
TTiF = +1 spin state distribution. The width of the images is 
1 mm. In f) the width of the mp = —1 component after a 
TOF = 25 ms, determined by a Gaussian fit, is plotted as a 
function of hold time (triangles). The fit to a Sigmoid func- 
tion (line) is included to guide the eye. The appearance of a 
small TTiF = il domains, which grow outward with hold time, 
contradicts a homogenous Bogoliubov theory, which predicts 
a uniform instability rate across the cloud. 

(see for Example IDd)). Since the mp = ^1 and mi? = 
clouds are immiscible for C2 > 0, this creates a poten- 
tial well that traps the pairs but allows the domain to 
grow axially due to the continued effect of the instability 
(expansion along the radial direction costs a substantial 
kinetic energy due to the tighter confinement). Thus the 
long timescale evolution of the instability exhibits one- 
dimensional coarsening dynamics We also observed 
smaller domain structures which could not be quantified 
clearly due to the presence of undamped mi? = density 
fluctuations in the initial state caused by nonadiabaticity 
in the initial transfer to the optical trap. 

In conclusion, we have tuned an antiferromagnetic con- 
densate through a phase boundary and quantified in de- 
tail the rate of instability in its vicinity. Future work 
will explore the phase coherence between the dynami- 
cally created mp = ±1 spin components in relation to 
topological defect formation. 
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